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ABSTRACT 

We report observations with the MIPS instrument aboard the Spitzer Space Telescope (SST) of four supernova 
remnants (SNRs) beUeved to be the result of core-collapse SNe: N132D (0525-69.6), N49B (0525-66.0), N23 
(0506-68.0), and 0453-68.5. All four of these SNRs were detected in whole at 24 fim and in part at 70 
fim. Comparisons with Chandra broadband X-ray images show an association of infrared (IR) emission with 
the blast wave. We attribute the observed IR emission to dust that has been collisionally heated by electrons 
and ions in the hot, X-ray emitting plasma, with grain size distributions appropriate for the LMC and the 
destruction of small grains via sputtering by ions. As with our earlier analysis of Type la SNRs, models can 
reproduce observed 70/24 pm flux ratios only if effects from sputtering are included, destroying small grains. 
We calculate the mass of dust swept up by the blast wave in these remnants, and we derive a dust-to-gas mass 
ratio of several times less than the often assumed value of 0.25% for the LMC. We believe that one explanation 
for this discrepancy could be porous (fluffy) dust grains. 

Subject headings: interstellar medium: dust — supernova remnants — Magellanic Clouds 



L INTRODUCTION 

Dust plays an important role in all stages of galaxy evolu- 
tion. The life-cycle of dust grains and the amount and relative 
abundances present in the interstellar medium (ISM) are de- 
termined by the balance betwe en dust form ation, grain modi- 
fication, and dust destruction ( Drainel200 3V Dust destruction 
is known to occur in both fast and slow shocks in SNRs ( Jon^ 
|2004). We focus here on dust destruction via sputtering by 
high energy ions in fast (non-radiative) shocks in SNRs. 

SNRs make excellent probes of the dust content of the dif- 
fuse ISM in galaxies, since their shock waves create X-ray 
plasmas that heat dust in their vicinities. Ambient dust heated 
by starlight may escape detection, but shock-heated dust ra- 
diates strongly in Spitzer wavelengths, so SNR studies probe 
different conditions than studies of diffuse ISM dust emission. 
Modeling of X-ray emission provides the basis for under- 
standing dust emission. Inferences of dust content and prop- 
erties from SNR observations are thus complementary to UV- 
absorption studies, which rely on the fortuitous locations of 
background UV-bright stars (Jenkins et al. 1984). The combi- 
nation of these types of investigation may lead to significant 
advances in our understanding of dust properties, with po- 
tential repercussions for gas-phase abundance determinations, 
theories of chemical evolution of galaxies, and dust-catalyzed 
cosmochemistry. 

To examine the nature of dust heating and destruction in 
the ISM, we conducted an imaging survey with the SST of 
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39 SNRs in the Magellanic Clouds. In a previous paper 
(Borkowski et al. 2006; Paper I) we analyzed IR emission 
from SNRs resulting from Type la SNe. We found that the 
dust-to-gas mass ratio was lower than the expected value 
of 0.25% (Weingartner & Draine 2001, hereafter WD), for 
the surrounding ISM, a result that we attributed to Type la 
SNe ex ploding in lower density media. However, Tap pe et alJ 
found a low dust content in N132D, a remnant of a 
core-collapse SN which occured in a dense ISM area, and sug- 
gested a high dust destruction efficiency in this SNR. A lower 
than expected dust content was also found in the circumstel- 
lar medium surrounding SN 1987A (Bouchet et al. 2006). In 
order to investigate this apparent dust deficit, we present here 
an analysis of four core-collapse SNRs, including N132D. 

We base our inference of the core-collapse nature of our 
objects on the presence of a pulsar-wind nebula in 0453- 
68.5 (Gaensleret al. 2003 ), a central co mpact object in N23 
dHavato et al.. .2006: Hushes et al."2006^, the O-rich classifi- 
cation of N132D (Lasker 1978), and the presence of a large 
mass of Mg in N49B (Park et al. 2003|). 

All four of these remnants were detected in both the 24 
and 70 /im bands of the Multiband Imaging Photometer for 
Spitzer (MIPS). We find a similar result to that found for the 
Type la remnants in Paper I: the dust-to-gas ratios we infer are 
lower bya factor of ^ 4 than what is expected for the LMC. In 
section 0] we discuss possible reasons for this apparent dust 
deficiency. 

2. OBSERVATIONS AND DATA REDUCTION 

All four objects were observed with the 24 and 70 pm MIPS 
arrays. At 24 /im, we mapped each remnant in our survey 
and the surrounding background with 14 frames of 30.93 sec- 
onds each, for a total exposure time of 433 s. At 70 /im, 
we observed a total of 545 s in 52 frames for all but N132D, 
which we observed for 986 s in 94 frames. MIPS images 
were processed from Basic Calibrated Data (BCD) to Post- 
BCD (PBCD) at the Spitzer Science Center (SSC) by ver- 
sion 13.2 of the PBCD pipeline. For the 70 /im images, we 
used the contributed software package GeRT to reprocess the 
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raw telescope images into BCD images, then reprocessed the 
BCD images using the SSC software package MOPEX. GeRT 
was useful in removing some of the artifacts, such as vertical 
stripes, from the 70 fim data. 

Figure 1 shows the 24 and 70 /im images as well as X- 
ray images from Chandra archival data and optical images 
from the Magellanic Cloud Emission-Line Survey (MCELS; 
Smith et al. 2005). Because line emission from low-ionization 
gases should be a significant contributor only in slower, radia- 
tive shocks, we conclude that we are seeing thermal IR emis- 
sion primarily resulting from dust. Recent spectroscopic ob- 
servations of N132D with Spitzer {Twppe et al. 2006) confirm 
that line contributions at 24 /im are negligible. Our measured 
fluxes are presented in Table[2 

3. MODELING 

Because of the morphological similarities between the IR 
emission and the blast wave seen in X-rays, we believe 
the dust present in the ISM is being collisionally heated by 
electrons and protons in the outward moving shock wave 
(|pwek & Arendt 1992). The modeling of dust emission for 
these remnants is similar to the modeling done on the four 
Type la remnants in Paper I, where it is described more fully. 
Our model uses a one-dimensional plane-shock approxima- 
tion. The sputtering timescale, Tp, is equal to J^ripdt, and 
is one of the inputs to the code. The other inputs are elec- 
tron temperature T^, ion temperature 7], gas density n, grain 
size distribution, and grain composition and relative abun- 
dances. For the composition, abundances, and distribution 
of dust grains, we follow WD, particularly the model con- 
sisting of separate silicate and carbonaceous grain popula- 
tions in fixed proportions, in a range of sizes from 1 nm to 
1 fim. We use the "average" LMC model with the maximum 
amount of small carbonaceous grains. We do not model emis- 
sion features from polycyclic aromatic hydrocarbons (PAHs), 
as these emission features are not expected to contribute in 
MIPS bands. In order to test the dependency of our results on 
the dust grain-size distribution, we also consider an alterna- 
tive model put forth by Cartledge et al. (2005), hereafter C05. 
Dust masses derived from this model were nearly identical to 
those from the WD model, but since the C05 model predicts 
a higher dust-to-gas ratio in the LMC (0.45%), this model ac- 
tually exaggerates the observed dust deficit. We thus report 
only the results from the WD model, and believe that these 
results are conservative. X-ray analysis provides estimates 
of the electron temperature, ionization timescale r,- = ngdt, 
and emission measure EM = UeMg (Mg is mass of swept-up 
gas) of the plasma. We assume tight dust-gas coupling, since 
the gyro-radii of dust grains are sma ll compared to the thick- 
ness of the shock jDwek et alHT99^ . We used archival Chan- 
dra data for our analysis, and fit X-ray spectra using Sedov 
non-eq uilibrium ionization (NE I) thermal models in XSPEC 
iArnaud. 1996,:.Borkowski et aL20 01 ). Emission-measure av- 
eraged Te and Ti from these models and the reduced r,- (1/3 of 
the Sedov Ti^ed, defined as the product of the postshock elec- 
tron density and the SNR age) were then used as inputs to our 
plane shock model. (The approximate factor 1/3 arises from 
applying results of a spherical model to plane-shock calcula- 
tions; see Fig. 4 of Borkowski et al. 2001). Since the Ti/Te 
ratio is close to 1 in these models, we set Ti = Te. 

We employed two methods to model dust emission. In the 
first, we fix Tp, Te, and 7} as derived from X-rays (taking 
Tp = Ti/ 1.2), leaving only the density and total dust mass as 
free parameters. In the second, we use the dynamical age of 



the remnant, as derived from optical or global X-ray studies, 
leaving the shock age and density as free, but correlated, pa- 
rameters. We find that these two different methods produce 
very similar results, within a factor of 25%. We thus report 
only the inputs and results from the first method. Model in- 
put parameters are given in Table|2l The density of the gas is 
then adjusted to reproduce the observed 70/24 /im flux ratio, 
and 24 fim flux is normalized to the observed value to pro- 
vide a total dust mass in the region of interest. The emission 
measure divided by the electron density gives an estimate for 
the amount of gas swept up by the blast wave, and dividing the 
dust mass by the gas mass gives us a dust-to-gas mass ratio for 
the shocked ISM. Finally, dust temperatures are calculated in 
our model. Because we use a range of grain sizes, we report 
a range of temperatures, since grains of different sizes will 
be heated differently under the same plasma conditions. The 
range of temperatures quoted is therefore the result of a single 
fit to the flux ratio. Results are summarized in Table|3] 

3.1. N132D 

SNRN13 2D has been well studied in optical wavelengths 
(Morse et al. 1996), and is one of the brightest remnants in 
the Magellanic Clouds at X-ray wavelengths. The remnant is 
extraordinarily bright at 24 nm, where we find a total flux of 

3 Jy, consistent with'Tappe e t al.N2006 t). It is the only rem- 
nant in our sample that is brighter at 24 /im than at 70 /im. 
This implies warm dust, and therefore a dense environment 
(to provide the inferred heating), consistent with the remnant 
being very bright in X-rays. Mor se et al.l {1596) estimate a 
preshock hydrogen density for this remnant of 3 cm"^ based 
on modeling of the photoionized shock precursor We ana- 
lyzed the NW rim and the bright southern rim separately. We 
find high densities, in good agreement with the expected post- 
shock proton density Hp of 12 cm"^. Densities are higher by a 
factor of ^ 2 in the NW, and comparatively less mass in gas 
in that region. The mass in dust, however, was comparable to 
what was found in the south, adjusting for the different sizes 
of the regions. 

3.2. N49B 

N49B is believed to be older than N132D, perhaps as old 
as 10,000 yrs. The ionization age is lower than this, how- 
ever, perhaps du e to a n explosion into a pre-existing cav- 
ity (Hushes et al.l ll998h . This remnant is much fainter than 
N132D at all wavelengths, but especially so at 24 /im, and the 
flux ratio is lower by more than order of magnitude in N49B 
compared to N 1 32D. This can be explained by the much lower 
density we find in N49B: = 1.1 cm"^ from our dust models, 
and Hp = 2.1 cm~^ from a Sedov analysis of X-ray data. Since 
dust heating rates are strongly dependent on density, a low- 
density environment will result in much cooler dust, whose 
spectrum peaks at longer wavelengths. Despite the contrasts 
in density and flux ratio of more than an order of magnitude, 
the dust-to-gas ratio in this remnant is only two times lower 
than that found in N132D. 

3.3. N23 

At 70 /im, the brightest parts of the shell of N23 are clearly 
visible. The region selected for analysis completely enclosed 
the emission visible at 70 /im, which contained about 75% 
of the emission visible at 24 /xm. Dust modeling yields rtp = 
5.8 cm""*. Hughes et al. (2006) found densities of 10 and 23 
cm"^ in a couple of X-ray bright shocks in this region of the 
remnant. 
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3.4. 0453-68.5 

Although 0453-68.5 had the highest 70/24 fj,m ratio of any 
remnant in this sample, only the north rim of the shell was 
clearly separable from the background at 70 fim. This is due 
to the overall faintness of the remnant, the high levels of back- 
ground in the region, and the lower sensitivity of the 70 fim 
array. Nonetheless, there was sufficient S/N to analyze the 
north rim of the remnant. The density found from modeling 
dust emission. Hp = 0.63 cm"-', was also the lowest of these 
remnants. Sedov modeling of X-rays yielded a value of Hj, 
= 1.1 cm"^. Both the temperature and density are consistent 
with the observations of overall faintness at 24 /xm and a high 
70/24 ratio. 

4. DISCUSSION AND CONCLUSIONS 

We find excellent (generally within a factor of 2) agreement 
between densities derived from dust modeling and estimated 

from opti cal and X-ray da ta. As has been proposed in the 

past (e.g., iDwek & A rendt 1992). dust emission is indeed a 
valuable density diagnostic of X-ray emitting plasmas. The 
good match between X-ray and IR derived densities suggests 
that the X-ray emitting gas is not highly clumped. 

Our basic quantitative results are contained in Table |3l 
First, note that the total mass in dust we infer is quite small, 
of order 0.01-0.1 Mq, for all objects. Morphological ev- 
idence suggests that this dust is associated with the blast 
wave, but these values also serve as restrictive upper lim- 
its for the amount of ejecta dust that can have been pro- 
duced in these core-collapse supernovae. They are, for the 
most part, less than the O.O8M0 < < 0.3Mq calculated 
bv llbdini & Ferrara (2001). Substantial amounts of dust 
can be destroyed when the ejecta are reverse-shocked, but 
the ejecta dust mass ultimately delivered to the ISM appears 
to be quite small. This is consistent with the findings of 
Isitanimirovic et al. (2005), who analyzed IR emission from 
the LMC SNR IE 0102.2-7219. 

We find that on average, about 40% of the mass in dust 
grains, including ^ 90% of the mass in grains smaller than 
0.04 /im, has been destroyed via sputtering in these remnants. 
After accounting for this sputtering, we find an initial dust-to- 
gas ratio that is lower by a factor of roughly four than what 
is generally expected for the LMC. The values in Table|3]are 
similar to the dust-to-gas ratios we found in remnants from 
Type la SNe, and they are still well below the value reported 
in WD of ~ 2.5 X 10"^. In Paper I we speculated that low 
dust/gas ratios might be expected for Type la remnants ex- 
panding into lower-than-average density material, in which 
previous supernovae might have destroyed some grains. Since 
we now have a similar deficit for core-collapse remnants ex- 
panding into denser media, this explanation is less likely, al- 
though it is still possible that the dust content is generally low 
in the vicinity of core-collapse SNe. 

The overall lower dust/gas ratios we infer will require an- 
other explanation. It is possible that sputtering timescales or 
rates might have been underestimated. Our models indicate 
that sputtering timescales would have to be increased by an 
order of magnitude, on average, to match the expected dust- 
to-gas ratio in the WD model, an unlikely possibility. Var- 
ious factors could increase sputtering rates. Sputtering can 
be enhanced by increased relative velocities between ions and 
grains, either caused by grain motions (not included in our 
models) or by higher than assumed ion temperatures. In or- 
der to estimate the importance of such effects, we doubled ion 



temperatures in our shock models, resulting in only 10-15% 
increase in the fraction of the original dust destroyed for the 
relatively fast shocks considered here. Sputtering is expected 
to be enhanced for nonspherical grains with increased surface- 
to-volume ratio S/V. For a prolate ellipsoidal grain with an 
axial ratio of 2, S/V is twice as large as for a spherical grain, 
which enhances sputtering rates by a modest factor of 2. We 
find it unlikely that these modest enhancements in sputtering 
rates can account for the lower than expected dust-to-gas mass 
ratios. 

Porous (fluffy) grains provide a possible explanation for 
the low dust-to-gas ratios. Porous grains (Mathis 1996) can 
be characterized by a porosity (fraction of grain volume de- 
void of material) V. Such grains have an extinction per 
unit mass generally larger than for spherical compact grains, 
f^shchinnikov et al. (2006) modeled visual and UV absorp- 
tion towards C, Oph and a Sco with a mixture of highly porous 
{P > 0.9) and compact grains with dust mass of only 70% and 
44% of the dust mass in the WD model. 

We expect enhanced sputtering rates for porous grains com- 
pared with spherical compact grains. A highly porous (V = 
0.9) grain has surface area 4.6 times larger than a compact 
grain of the same mass, so for the same sputtering yields an 
enhancement of sputtering rates by a factor of 4.6 is expected. 
In addition, protons, a particles, and heavier ions can pen- 
etrate much deeper into porous grains, because their range 
within the solid is inversely proportional to the mean grain 
density. A proton with energy E will penetrate 0.09(£' / 1 keV) 
fim into a porous silicate grain with V = 0.9 (we used eq. 18 
of Draine & Salpeter 1979 to estimate the proton projected 
range Rp). Jurac et al. ( 1998) found that sputtering yields are 
enhanced if the grain radius a is less than 3Rp because of addi- 
tional sputtering from the back and sides of the grain. Grains 
as large as 0.25(£'/l keV) fim are then destroyed more ef- 
ficiently. It is also possible that the rough surface geometry 
expected for porous grains could result in a further enhance- 
ment of sputtering rates, as the protruding grain extensions 
comparable in size to the proton or a particle projected ranges 
could be more readily sputtered dJurac et alj ri998). An order 
of magnitude enhancement in sputtering rates is thus possible 
for highly porous grains. 

The reduced preshock dust content and substantial enhance- 
ment in sputtering rates that result from porous grain models 
provide a promising explanation for the apparent dust deficit 
in LMC SNRs. These effects might be partially offset by the 
increa sed IR emissivity of porous grains ( Voshchinnikov et aQ 
as less dust is then required to account for the ob- 
served IR emission. It is possible, however, that compact 
grains dominate the observed IR emission, if porous grains (or 
porous mantles surrounding compact cores) have been pref- 
erentially destroyed. Modeling of Spitzer observations with 
porous grain models is required to assess whether these mod- 
els can resolve the apparent dust deficit in LMC SNRs, includ- 
ing SNR 1987 A. Spectroscopic follow-up on these remnants 
in the far-IR is required to confirm and strengthen our results. 
Implications of a resolution go beyond these particular SNRs 
or even the LMC as a whole and could help bring about a sig- 
nificant improvement in our understanding of interstellar dust 
in general. 

This work was supported by NASA through Spitzer Guest 
Observer grant RSA 1265236. We thank an anonymous ref- 
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TABLE 1 
Measured Fluxes 



Object 24 /xm 70 fim 70/24 

N132DNW 730 ±73 430 ±96 0.59 ±0.13 

N132S 1000 ±100 770 ±170 0.76 ±0.17 

N49B 43 ±4.3 395 ±79 9.1 ±2.0 

0453-68.5 13 ±1.3 250 ±50 19 ±4.2 

N23 100 ±10 240 ±50 2.4 ±0.5 



"All fluxes given in millijanskys 



TABLE 2 
Model Inputs 



Object 


Tei= T.) 




Age 




(keV) 


(10^ cm"3 yr) 


(yr) 


N132D NW 


0.61 


17 


2500" 


N132DS 


1.02 


3.2 


2500 


N49B 


0.36 


2.4 


10,900'' 


0453-68.5 


0.29 


8.6 


8700'' 


N23 


0.56 


6.7 


4600*^ 



';Morse et al. ( 1996) 
' Gaensler et al j t2006i^ 
•[Hughes et alj 120061) 



TABLE 3 
Model Results 



Object 


lip (cm 


T(dust) (K) 


Dust Mass (Mq) 


iieMg (cm ^Mq) 


Dust/gas 


% dest. 


Dust/Gas (orig.) 


Lib 


N132D NW 


34 


95-120 


0.0075 


660 


4.6 X 10"^ 


50 


9.2 X 10-^ 


35 


N132D S 


14 


85-105 


0.015 


1020 


2.5 X 10"^ 


38 


4.0 X 10^ 


94 


N49B 


1.1 


50-60 


0.08 


580 


1.8 X 10"^ 


27 


2.4 X 10-4 


8.4 


0453-68.5 


0.63 


40-55 


0.10 


120 


6.1 X 10-^ 


33 


9.8 X 10~* 


4.5 


N23 


5.8 


70-85 


0.011 


250 


3.0 X 10"^ 


39 


4.9 X 10-4 


8.7 



Note. — Col (2): postshock proton density; col. (3): for 0.02-0.1 /im grains; col (4); mass of dust currently observed, after sputtering (post-shock), in Mq; 
col (6): ratio of current dust mass to swept-up gas mass; col (7): percentage of dust destroyed via sputtering; col (8): inferred dust-to-gas mass ratio in preshock 
ISM; col (9): L36 EE L,r/10^^ erg s"'. 




Fig. 1 . — Top row, from left to right: N132D at 24 and 70 /im (the region of interest at 70 microns is marked on the image), in the X-rays (broadband, Chandra 
image) and in the optical (overlay of MCELS images, with [S II], Ha and [O III] marked in red, green and blue, respectively). Second, third and fourth rows 
show the same sequence for N49B, N23 and 0453-68.5, respectively. One arc-minute scales are shown on the 24 fim images. 



